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ABSTRACT 

At redshifts beyond z > 6, as the mean fraction of neutral hydrogen (xjjj) in the intergalactic medium 
(IGM) increases, the line flux of Lya emitters can be significantly suppressed, which can result in a 
decrease in the observed number of emitters at a given Lya flux. However, cosmological HII regions 
surrounding the Lya emitting galaxies alleviate these effects. We use simple models of the Lya line sup- 
pression that incorporate the presence of HII regions to predict the overall effect of the Lya absorption 
on the Lya luminosity function. We find, in agreement with other recent studies, that unless ionizing 
sources are unusually strongly clustered, a fully neutral IGM may be inconsistent with the large abun- 
dance of confirmed z = 6.5 Lya emitters. However, the presence of local HII regions prohibits placing a 
tight constraint on the mean neutral fraction. We find (xhi) 0.25; the presence of strong winds and/or 
the clustering of ionizing sources would further weaken this constraint. We conclude that the evolution 
of the Lya LF is consistent with reionization occurring near this redshift. Finally, we suggest that a 
measurement of observed Lya line width as a function of the Lya luminosity, in a future, larger sample 
of Lya emitters, may serve as a robust diagnostic of the neutral fraction in the IGM. 

Subject headings: cosmology: theory - cosmology: early universe - galaxies: formation - galaxies: high- 
redshift - galaxies: quasars: general - galaxies: quasars: absorption lines 



Of 

6 

CO 



X 



1. INTRODUCTION 

The recent detection of Gunn-Peterson (GP) troughs 
in the spectra of the SDSS quasars (Becker et al. 2001; 
White et al. 2003) places a relatively direct lower limit on 
the mean (volume averaged) neutral fraction of the IGM 
at z ~ 6. While this limit is relatively weak, (xhi) 25 10~ 3 , 
the spectral imprint (Mesinger & Haiman 2004) and size 
(Wyithe & Loeb 2004a) of cosmological HII regions around 
the SDSS quasars indicate that the neutral fraction around 
z ~ 6 is significantly higher, (xhi) 0.1. The rapid evo- 
lution of (xhi) with redshift, inferred from the spectra of 
several quasars between 5.5 ^ z <S 6.5 (Fan et al. 2002; 
Cen & McDonald 2002; Lidz et al. 2002; but see Songaila 
2004), and also the thermal state of the IGM at 3 <J z <, 4 
(Hui & Haiman 2003), indicate and that a percolation of 
discrete HII regions is occurring around this epoch. 

Lya emission lines from high-redshift sources can serve 
as another probe of the ionization state of the IGM. The 
damping wing of the GP absorption from the IGM can 
cause a characteristic absorption feature (Miralda-Escude 
1998). In a significantly neutral IGM, the absorption can 
produce conspicuous effects, i.e. attenuating the emission 
line, making it asymmetric, and shifting its apparent peak 
to longer wavelengths (Haiman 2002; Santos 2004). These 
effects generally increase with (xhi), and it has been sug- 
gested that a strong drop in the abundance of Lya emitters 
beyond the reionization redshift can serve as a diagnostic 
of a nearly neutral IGM (Haiman & Spaans 1999). Re- 
cently, Malhotra & Rhoads (2004; hereafter MR04) and 



Stern et al. (2004) carried out the first application of this 
technique, by comparing the luminosity functions (LFs) of 
Lya emitters at z = 5.7 and z = 6.5. The LF shows no 
evolution in this range, and this has been interpreted as 
evidence against percolation taking place near z ~ 6. 

A potential caveat for this interpretation is that z ~ 6 
galaxies will be surrounded by their own local cosmologi- 
cal HII regions, which can significantly reduce the atten- 
uation of the Lya line flux (Cen & Haiman 2000; Madau 
& Rees 2000). It has been shown that this can render 
Lya lines detectable even if the galaxies are embedded in 
a fully neutral IGM (Haiman 2002; Santos 2004; Cen et 
al. 2004), especially when the increased transmission due 
to the clustering of ionizing sources is taken into account 
(Furlanetto et al. 2004; Gnedin & Prada 2004; Wyithe & 
Loeb 2004b), and/or if the Lya emission has a significant 
recession velocity with respect to the absorbing gas. 

In this Letter, we study the impact of IGM absorption 
on the Lya LF. Our analysis is similar to that of the re- 
cent analysis by MR04; the improvement is that we model 
the attenuation of Lya lines including cosmological HII 
regions. While very little attenuation (by a factor of <, 2) 
can be tolerated by the lack of evolution of the LF over the 
range 5.7 z 55 6.5, we find that the presence of local HII 
regions allows a neutral fraction as high as (xhi) = 0.25. 
On the other hand, (xhi) ~ 1 is allowed only if the ion- 
izing sources are unusually strongly clustered. In either 
case, we argue that the present Lya LFs are consistent 
with reionization occurring near z ~ 6. 1 



throughout this paper we assume a cold-dark matter cosmology (ACDM), with (Oa, ^6, ^o) = (0.73, 0.27, 0.044, 71 km s 1 Mpc 1 ), 
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2. THE OBSERVED LYa LUMINOSITY FUNCTION 

Following Malhotra & Rhoads (2004; hereafter MR04) 
we take as basic inputs to our analysis the Lya LFs at 
z = 5.7 and z — 6.5. We adopt the forms of the LF at 
these two redshifts from MR04. Note that at both red- 
shifts, the LFs are inferred by culling several independent 
datasets (Hu et al. 2002; Rhoads et al. 2003; Kodaira 
et al. 2004). At z = 5.7, the LF is comparatively well 
determined and well fit by a Schechter function with rela- 
tively small errors in normalization (<E>+) and characteris- 
tic luminosity (L*), while the faint-end slope (a) is less 
well-determined. Here we adopt the best-fit values of 
($*, -£,*,«) = (10~ 4 Mpc- 3 ,10 43 erg s-\-1.5) (displayed 
by the short-dashed curve in Fig. 2 below). Within the 
uncertainties quoted in MR04, the values of these param- 
eters do not significantly change our conclusions below. 
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Fig. 1. — The suppression of the total Lya line flux, as a 
function of the observed line luminosity, for an IGM with 
(xhi) = 0.1 (or Ti2 = 0.0015; long-dashed curve), and for 
and IGM with (x H i) = 1-0 (or T 12 = 0; solid curve). 

At z = 6.5, the LF is less well determined, as shown 
by four data-points with errors in Figure 2 below. Of 
these four data-points, the most constraining is from the 
deep Subaru field, which revealed the presence of a hand- 
ful of emitters at z ~ 6.5 (Kodaira et al. 2003), and has 
the smallest error bar. The single, faint lcnsed galaxy at 
this rcdshift found by Hu et al. (2002; represented by the 
lowest-luminosity point in Fig. 2 below) implies a surpris- 
ingly high abundance of the faintest emitters, but with a 
single source, it has a large uncertainty. Finally, we fol- 
low MR04 and conservatively assume that the intrinsic LF 
(i.e. prior to processing of the line through the IGM) is 
the same at z = 6.5 as at z — 5.7. More realistic models, 
based on the hierarchical growth of structures, combined 
with Lya line processing by gas and dust interior to the 
galaxies, predict that the number of emitters were smaller 
at earlier epochs (Haiman & Spaans 1999; Le Delliou et al. 



2004), which would strengthen the limits on (xhi) below. 

3. THE ATTENUATION OF Lya EMISSION LINES 

We follow Cen & Haiman (2000), and use a simple model 
to find the attenuation of the Lya line as a function of 
wavelength. The model is straightforward, and we only 
briefly recap the main features. We start with a Gaussian 
emission line originating at z s — 6.5, with a line-width of 
Aw = 300 km s _1 , typical of normal galaxies (narrower 
lines would be attenuated more severely; see discussion 
below). We assume that the source is surrounded by a 
spherical Stromgren sphere that propagated into the IGM 
with a neutral fraction (xhi)- The IGM is assumed to 
have a clumping factor C = (p 2 )/(p) 2 = 10 with a log- 
normal distribution. We assume the presence of a uniform 
ionizing background, characterized by an ionizing rate Ti2 
(ionizations per 10 12 s per atom, yielding the given value 
of (xhi) in ionization equilibrium for the IGM. 

We assume a Salpeter IMF, and relation (Kennicutt 
1983) between star-formation rate and Lya luminosity 
(SFR=1 M Q /yr) (L(Lya) = 10 42 erg/s). We conser- 
vatively assume / osc = 1 (escape fraction of ionizing ra- 
diation) in computing R s , and a lifetime t s — 10 8 yrs for 
each source. With all other parameters in the model fixed, 
we vary Fi 2 (or equivalently (xhi)), and compute the fac- 
tor by which the total line flux is suppressed relative to 
the unabsorbcd Gaussian line. This factor is shown in 
Figure 1 as a function of the intrinsic luminosity of the 
source. Faint sources are more severely attenuated: by a 
factor of up to ~ 15 for log L (Lya /erg s _1 ) =40 [or up to 
~ 6 for logi(Lya /erg s _1 ) = 42, the present lower limit 
for the detection of a Lya emitter at z = 6.5] in a neutral 
universe. Bright (logL(Lya /erg s _1 ) = 43) sources, how- 
ever, are only attenuated by a factor of ~ 1.5 even if (xhi) 
is as high as 0.1. 

The attenuation factors in Figure 1 are computed by as- 
suming that the HII regions are being created by a single 
ionizing source (the Lya emitter itself). These HII regions 
(with typical sizes of a few Mpc in the observed range of 
source luminosities) may contain many other undetected 
galaxies, which make the HII regions larger (Furlanetto et 
al. 2004; Gnedin & Prada 2004; Wyithe & Loeb 2004b) 
and also more highly ionized. Therefore, the results shown 
in Figure 1 should be viewed as conservative upper limits. 
These results show that the HII regions of individual galax- 
ies can significantly reduce the attenuation, especially for 
relatively bright sources. 

4. THE IMPACT ON THE LYq LUMINOSITY FUNCTION 

We next use the attenuation factors computed above, 
and obtain the reduction in the cumulative Lya LF. The 
resulting LF is shown in Figure 2, for three values of (xhi): 
0, 0.1, and 1 (top to bottom). The plot also shows the 
data for the z = 6.5 LF with error bars (from MR04). 
The short-dashed curve is a reproduction of the z — 5.7 
LF (a Schechter function), which, by assumption, would 
be identical to the z = 6.5 LF if the IGM contained no 
neutral HI. 

Note that the suppression of the line flux depends on 
the luminosity of the source, but the reduction of the LF 
in Figure 2 is uniform, by a nearly constant factor. The 



consistent with the recent results from WMAP (Spergcl et al. 2003). Unless stated otherwise, all lengths are quoted in comoving units. 
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steepening of the LFs at the bright end compensates for 
the reduced effect of line attenuation; the two conspire to 
give a nearly constant suppression of the cumulative LF. 
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Fig. 2. — The Lya luminosity functions for a fully ionized 
IGM (no suppression, short-dashed curve), for an IGM 
with (xm) = 0.1 (or T12 = 0.0015; long-dashed curve), 
and for and IGM with (xm) = 1-0 (or Ti 2 = 0; solid 
curve) . 



5. RESULTS AND DISCUSSION 

We can now assign a likelihood to each value of (xm), by 
comparing the model LFs and the data shown in Figure 2 . 
For each value of (xm), we obtain the mean expected num- 
ber of Lya emitters above the flux thresholds shown by the 
data points (here we use the effective volume probed by 
each survey, listed in Table 2 in MR04). We assume Pois- 
son fluctuations to compute the likelihood in each bin. We 
find that the lowest-luminosity data-point has an a-priori 
Poisson likelihood that is exceedingly low, even in a fully 
ionized universe. As Figure 2 shows, <, 0.03 objects are 
expected in the small survey volume in which the single 
lensed Lya emitter, represented by this data-point, was 
discovered. However, given the difficulties in obtaining an 
accurate effective volume from the poorly constraint lcns- 
ing configuration for this object, we follow MR04 and omit 
this source from our analysis. The result, i.e. the product 
of the Poisson probabilities for the three higher luminosity 
bins, are shown in Figure 3 as a function of (xm)- 

We find that (xm) = 1 is ruled out at high significance, 
in agreement with the conclusions reached by MR04 and 
Stern ct al. (2004) (but see discussion below). On the 
other hand, at the 99.7% confidence level, we find only 
the relatively weak constraint (xm) < 0.25. The upper 
dashed line corresponds to the limit of no HI absorption 
and shows that the data is consistent (at the ~ 10% con- 
fidence level) with the intrinsic z = 5.7 and z — 6.5 LFs 
being identical (but note that the IGM would have to be 
ionized to an unrcalistically high level to actually approach 
this limit). 
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Fig. 3. — The probability of each mean neutral fraction 
in the IGM, given the observational constraints from the 
z = 6.5 Lya luminosity function. The model assumes that 
the intrinsic LF is the same as at z — 5.7 (see text for 
discussion). 

As mentioned above, HII regions can be enlarged due 
to the clustering of ionizing sources, and also due to their 
being elongated along our line of sight (e.g. Gnedin & 
Prada 2004). These effects reduce the Lya attenuation, 
and weaken the upper limit on (xhi)- It is interesting to 
ask whether any constraint can be placed on the neutral 
fraction, in the presence of clustering. To address this 
question, we re-computed the attenuation factors shown 
in Figure 1, and the implied LFs shown in Figure 2, fix- 
ing (xhi) = 1, but boosting the ionizing emissivity of each 
source by a constant factor of B. We find that the boost- 
ing factor required to render this model acceptable at the 
99.7% confidence level is approximately B = 60 (note that 
the increased luminosity is required to reduce both the GP 
damping wing, and the residual HI inside the HII region). 

In order to assess whether the ionizing emissivity can be 
boosted by this factor, we computed the total mass in dark 
matter halos in the mass range M m i n to M s and within a 
sphere of radius R around a halo of mass M s . We included 
the spatial correlations between halos with a prescription 
for linear bias (see eq. 15 in Haiman ct al. 2001). For the 
likely halo mass range of 10 n M Q <, M <, 10 12 M Q (Rhoads 
et al. 2003) and HII sphere radii (2 <, R <, 6 Mpc; e.g. 
Haiman 2002) of Lya emitters at z = 6.5, we find that the 
total collapsed halo mass, down to 2 x 1O 8 M0 (correspond- 
ing to a virial temperature of 10 4 K at z = 6.5) is increased 
by a factor of 2 — 18 (the satellite halos have linear bias 
parameters of 2 <; b <^ 10). Unless the star formation 
efficiency is preferentially higher in the low-mass "satel- 
lite" galaxies, it is unlikely that the ionizing emissivity is 
boosted by the required factor of 60. A recent numerical 
simulation by Gnedin & Prada (2004, see their fig. 3) how- 
ever, suggests that at least 10% of galaxies may transmit 
their Lya lines without any absorption at z = 6.5. 
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Fig. 4.— The FWHM of the transmitted Lya line as 
a function of the luminosity of the source. The upper 
(lower) panel assumes an intrinsic Gaussian line shape 
with a width of 300 (30) km/s. On both panels, the 
dashed curve assumes an ionizing background of Ti 2 = 0.1 
(or (xm) ~ 10~ 3 ), and the solid curve assumes Ti 2 = 
((arm) = !)• The feature at - 10 41 ~ 42 
to a significantly neutral IGM. 
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Another effect that can further weaken the upper limits 
on the neutral fraction is a systematic redshift of the Lya 
emission lines relative to the absorbing gas. Lyman break 
galaxies at lower redshift reveal such systematic shift by 
several hundred km/s, attributed to galactic winds (Shap- 
ley et al. 2003). If velocity offsets as high as this are 
common for Lya emitters at z ~ 6, this would render Lya 
line attenuation from the residual HI within the HII region 
effectively negligible, and will also significantly reduce the 
GP damping wing absorption. By repeating our analysis 
above to obtain a likelihood for a neutral universe, we find 
that a redshift by 600 km/s is required to for (xhi) = 1 to 
yield an acceptable fit for the LF. Shapley et al. (2003) 
find a decreasing velocity offset with increasing Lya equiv- 



alent width, with {Av} < 500 km/s for W a > 50A (see 
their fig. 11), suggesting that high-z Lya emitters may 
not have the requisite 600 km/s offsets. 

Finally, we propose a different diagnostic of the neutral 
fraction that could be available from a future, larger sam- 
ple of z 6 Lya emitters. For faint Lya emitters, the line 
attenuation is dominated by the residual HI inside their 
cosmic HII regions, whereas for bright emitters, the GP 
damping wing is dominant (Haiman 2002). In Figure 4, we 
show the predicted full width at half maximum (FWHM) 
of the transmitted line, as a function of Lya luminosity 
(for two different intrinsic line-widths). The figure reveals 
a clear imprint of this transition at ~ 10 41 ~ 42 ergs -1 . 
This feature occurs only if the IGM contains a significant 
amount of neutral hydrogen (solid curves), and we checked 
that it is present for different assumed intrinsic line-shapes 
(e.g. a top-hat or a Lorentzian, which have less/more 
extended wings than a Gaussian). This suggests that a 
measurement of observed Lya line width as a function of 
the Lya luminosity may therefore serve as a more robust 
diagnostic of the neutral fraction. 

6. CONCLUSIONS 

In this Letter, we have considered the attenuation of the 
line flux of Lya emitting galaxies at z = 6.5. We included 
the effect of HII regions around these galaxies, which re- 
duce the line attenuation. We computed the overall impact 
of the IGM on the Lya luminosity function, and found that 
the presence of the HII regions makes a large mean neutral 
fraction in the IGM acceptable. The best upper limit is at 
most as strong as (xm) Si 0.25. However, a fully neutral 
IGM may be acceptable only in the presence of an unusu- 
ally strong clustering of ionizing sources, c md/or a strong 
systematic recession velocity of the Lya emitting gas. In 
either case, the present determination of the z ~ 6.5 Lya 
LF is consistent with reionization occurring at this epoch. 
The LF of a larger sample of Lya emitters, or a systematic 
study of their linewidths and shapes, as a function of lu- 
minosity, should nevertheless provide a robust diagnostic 
of the ionization state of the IGM in the future. 
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